Topological phase transitions in spin-S ladders 
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We show that spin-S ladders undergo, for antiferromagnetic rung coupling, a series of 4S topo- 
logical phase transitions when the ratio of the leg to rung coupling changes from strongly antiferro- 
magnetic to strongly ferromagnetic. These phase transitions are characterized by a change of Berry 
phase on the rungs, and, with appropriate boundary conditions, by the appearance and disappear- 
ance of edge states without the bulk gap ever closing. The various phases can be distinguished by 
the number of valence bond singlets on the rungs, which decreases from 2S for strong rung coupling 
to for dominant leg coupling. In particular, the isotropic spin-3/2 ladder with equal leg and rung 
coupling is explicitly shown to be in the valence bond solid phase with one singlet per bond. 

PACS numbers: 75.10.Jm,75.10.Pq,75.40.Mg 



Topological aspects of matter have become one of the 
dominant themes in solid state physics. The character- 
ization of quantum phases by topological invariants has 
been an extremely fruitful concept in the Quantum Hall 
Effect [U [2], and it lies at the root of more recent devel- 
opments such as topological insulators |3l |4]. A related 
question concerns the nature of the phase transitions be- 
tween quantum phases with different topological proper- 
ties. Quantum phase transitions [S] that involve a sym- 
metry breaking are often continuous, and, even if both 
phases are gapped, the gap closes at the transition. By 
contrast, transitions that do not involve any symmetry 
breaking are usually first-order, and it is in principle pos- 
sible to go continuously from one phase to the other by 
taking another path, as in the liquid-gas transition by 
going around the critical point. In the case of transitions 
between phases that are only distinguished by topologi- 
cal properties, much less is known. It has been recently 
argued that, for topological phases to be really differ- 
ent, the system must possess some symmetry to protect 
them[6|. If this is the case, i.e. if phases are really dis- 
tinct, there is however no general theory of what hap- 
pens at the transition. In particular, under which cir- 
cumstances should one expect edge states to appear and 
the bulk gap to close is not known. Simple examples 
of topological phase transitions are clearly called for to 
make progress in understanding their general properties. 

In this Letter, we show that spin-S ladders with anti- 
ferromagnetic rung coupling undergo a cascade of topo- 
logical phase transitions characterized by a topological 
invariant and by the appearance and disappearance of 
edge states in a parameter range where the bulk gap 
never closes, and where, accordingly, no quantum phase 
transition takes place. Spin ladders are defined by the 
Hamiltonian 



H 



Ladder 



— J\\/_^ /_^ Si_Q,-Si+i^ct + Jl ^^ Si^i •Si_2 (1) 



where i is the rung index, a the leg one, and S^^q are spin 
S operators. Throughout we will use the parametrization 



J|| = J cos 9 and J_l — J sin 6*. Over the last two decades, 
spin ladders have been investigated in great detail 1^. It 
has been shown early on that, for antiferromagnetic rung 
coupling, half-integer spin ladders are gapped if the num- 
ber of legs is even and gapless if it odd while integer spin 
ladders are always gapped, and for the two-leg ladder, 
a consensus has emerged, based on numerical investiga- 
tions and field theory arguments, that there is no quan- 
tum phase transition between the weak and strong rung 
coupling regimes. This is to be contrasted with the case 
of dimerized spin chains, in which a series of quantum 
phase transitions take place upon increasing the dimer- 
ization, but despite specific investigations aimed at find- 
ing similar quantum phase transitions in the spin-1 two- 
leg ladder [S], the evidence has imposed itself that the spin 
gap remains open all the way from weak to strong rung 
coupling, definitely excluding the presence of a quantum 
phase transition. The only phase transition that has been 
detected in spin-1/2 ladders with antiferromagnetic leg 
coupling takes place when the sign of the rung coupling 
changes from antiferromagnetic to ferromagnetic. At this 
transition, the gap closes, and the two topologically dis- 
tinct singlet phases can be distinguished by the type of 
string-order parameter (even or odd) that exhibits long- 
range order [5HTT]. 

The starting point of the present study is the inves- 
tigation of the Berry phase on the rungs of the lad- 
ders. The Berry phase can be defined for any Hamil- 
tonian H{<f)) which depends periodically on a parameter 
(/)fT^. If \GS{(f))) denotes a single-valued ground state 
of H {(})), the Berry connection is defined by A{(f)) — 
{GS{(t))\d(i,\GS{(j))), and the Berry phase is the integra- 
tion of the Berry connection over a loop: ij — ^ A((j))d(l). 
A few years ago, Hatsugai [T^ has shown that one can 
detect the presence of a singlet on a given bond of an 
antiferromagnet by calculating the Berry phase asso- 
ciated to a twist of the transverse component of the 
spin-spin interaction on this bond: S^S^ + S^ SJ' — > 
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Figure 1: (Color online) Phase diagrams of the spin-S ladders for S=l/2, 1, 3/2 and 2 as a function of 9 with the convention 
J± = J sin 6' and Jy = J cos 6. They consist of four types of phases: a ferromagnetic (FM) phase (tt < 6 < 3n/2), a Haldane 
phase (37r/2 < 9 < 2n), and two types of singlet phases for antiferromagnetic rung coupling (0 < 6^ < tt). These singlet phases 
are characterized by a topological invariant, the Berry phase on a rung, which is equal to in the light (yellow) phases and 
TT in the dark (blue) phases. The integer number reported in each phase corresponds to the number of valence-bond singlets 
according to the valence-bond interpretation (see main text). It takes values between to 2S, leading to a total of 4S phase 
transitions for antiferromagnetic rung coupling. For S=l/2, 1 and 3/2, the critical values of 9 are the result of a finite-size 
analysis, while for S=2 they correspond to the values for 8 sites. 
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In a valence bond solid, the Berry phase is related to 
the number Ng of singlets of the bond by: 



7 = ^^. 



(mod 27r), 



(2) 



this quantization being protected by time-reversal sym- 
metry. 

This criterion has been tested on dimerized spin chains 
with larger spins [Tl], which undergo a series of quantum 
phase transitions upon increasing the dimerization|15j. 
and indeed the Berry phase jumps between and tt at 
each transition, in agreement with the interpretation of 
the phases in terms of valence-bond singlets. 

For spin ladders, since there is no quantum phase 
transition upon increasing the rung coupling, one might 
naively expect that the Berry phase does not change upon 
increasing the rung coupling, and that it is always equal 
to the value it has to take in the limit of strong rung 
coupling, TT for half-integer spins with an odd number 
of valence-bond singlets and for integer spins with an 
even number of valence-bond singlets. To investigate 
this point further, we have calculated the Berry phase 
for spin-S ladders with S=l/2, 1, 3/2 and 2 as a func- 
tion of the rung coupling (see Supplemental Material for 
details). The calculations have been performed on fi- 
nite ladders with periodic boundary conditions (PBC). 
It turns out that, for antiferromagnetic rung coupling, 
the Berry phase changes several times between and 
TT upon increasing the rung coupling for both antiferro- 
magnetic and ferromagnetic leg coupling, resulting in the 
phase diagrams of Fig. [T] The critical values of 9 have 
been deduced from a finite scaling analysis, except for 
S=2, where results for 8 sites are shown. The scaling is 
quite smooth, up to an even-odd effect for antiferromag- 
netic leg coupling, as can be seen from the example of the 
spin-1/2 ladder (see Fig. [2]). The finite-size scalings for 
spin-1 and spin-3/2 ladders are shown in the Supplemen- 
tal Material. Note that the transition a.t 9 = 0.26 tt for 



spin 1/2 has already been reported in Ref. 16 in the con- 
text of an investigation of the frustrated spin-1/2 ladder 
with four-spin ring exchange. Note also that this tran- 
sition does not show up in the string order parameters 
that distinguish the phases with ferromagnetic and anti- 
ferromagnetic rung couplings [5Hl 1) . 

To be complete, we have also included the ferromag- 
netic rung coupling part of the phase diagram, which sim- 
ply consists of two phases: a ferromagnetic phase when 
the leg coupling is ferromagnetic (see Supplemental Ma- 
terial), and a Haldane phase [T7] when the leg coupling 
is antiferromagnetic. For strong ferromagnetic rung cou- 
pling, the spins of a rung bind into a spin 2S, and the 
low-energy effective model is just a spin-2S Heisenberg 
chain. This Haldane chain can be characterized by the 
Berry phase on one of its bonds, which is equal to n for 
2S odd and for 2S even. For the ladder, this requires 
to calculate the Berry phase associated to simultaneous 
twists on two leg bonds facing each other. We have found 
that this Berry phase is constant throughout that phase, 
with a value equal to tt for half-integer spin (2S odd) and 
to for integer spin (2S even), in agreement with the 
interpretation of this phase as a Haldane phase. 

For — 7r/2 < 9 < n, the spin gap never closes for in- 
teger spins, and it only closes at the transition between 
the Haldane phase and the other singlet phases for half- 
integer spins (for 9 = 0, the system consists of decoupled 
half- integer spin chains, which are gapless). So the phase 
transitions detected by the Berry phase for < 9 < tt 
cannot be quantum phase transitions. Since the Berry 
phase is a topological invariant, they must be topological 
in nature. The next question is thus to determine the 
nature of the various phases. In Ref. [16] , the authors 
proposed the name "rung * singlet" for the weak-rung 
phase below 9 = 0.26 tt, suggesting that it is another 
type of rung singlet phase to be further investigated. 
However, in view of the cascade of phase transitions for 
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Figure 2: (Color online) Finite-size scaling of the topological 
transitions of the spin-1/2 ladder. The results have been fitted 
with a fourth order polynomial in l/N . 



larger spins, it seems more natural, by analogy with the 
dimerized chains, to interpret the various phases as cor- 
responding to different numbers of valence-bond singlets 
on the rungs, going from 2S in the strong rung limit, 
when pairs of spins form singlets on each rung, to upon 
approaching the limit of decoupled chains. According to 
this interpretation, the weak-rung phase of the spin-1/2 
ladder is not a rung singlet phase since it has no singlet 
on the rung. 

To check this interpretation in terms of valence-bond 
singlets, we have investigated another possible character- 
istics of topological phase transitions, namely the appear- 
ance or disappearance of edge states at the transition. In 
the Haldane phase, it has long been known that edge 
states appear in finite chains [THHH]. In the present case, 
these edge states are expected to appear if we consider 
finite ladders with vertical edges, and indeed we have ob- 
tained numerical evidence of the presence of such edge 
states in the Haldane phase of the spin-1/2 ladder (not 
shown). By contrast, to distinguish the singlet phases 
for antiferromagnetic rung coupling, one should rather 
consider open boundary conditions (OBC) with diagonal 
edges [22], Indeed, in the strong rung limit, the spins on 
the rungs are bound into singlets, and the orphan spin 
at each end of the ladder is expected to give rise to edge 
states. So we have investigated the fate of the edge states 
with diagonal edges as a function of the rung coupling in 
a spin-1/2 ladder. The low energy states with PEC and 
with diagonal edges are plotted as the inverse number of 
spins for two values of theta: 9 = 0.3 tt, where the Berry 
phase is equal to 1, and 9 = 0.1 tt, where the Berry phase 
is equal to 0. For PBC, the gap is finite in all cases, in 
agreement with previous results. For OBC with diagonal 
edges, edge states are indeed present for 9 = 0.3 tt, but 
they disappear when the Berry phase changes from tt to 
0, and they are not present any more for 9 = 0.1 vr. This 
confirms the interpretation of the change of Berry phase 
as a topological phase transition, with the disappearance 
of edge states upon reducing the rung coupling. 

What is remarkable as compared to other cases where 
edge states appear is that there is neither a first-order 



transition (unlike the case of frustrated S* > 1 spin 
chains p3l I24j). nor a continuous phase transition with 
a gap closing, a standard argument to explain the ap- 
pearance of edge states at a topological phase transition 
relying on the gap closing at the transition upon going 
from one phase to the other. So, to better understand 
the fate of the edge states at the transition, we have cal- 
culated the evolution of the magnetization profile in the 
Sz = ^ sector of a ladder with diagonal boundary con- 
ditions from the strong rung limit to the transition at 
9 = 0.26 TT. In the strong rung limit, the magnetization 
is concentrated at the edges of the system, and the fast 
decay can be fitted with an exponential with a very small 
length scale ^*^ (see Supplemental Material). Upon ap- 
proaching the transition, the magnetization is less and 
less localized, a behavior consistent with a divergence of 
^^^. A similar information can be extracted from the en- 
ergy spectrum: the energy of the edge states collapse onto 
the ground state exponentially with the inverse number 
of rungs, and the length £^^ that can be extracted from 
the fit is in good agreement with ^^^. 

So, at the transition, a typical length scale seems to 
diverge. This is usually associated to a gap closing, but 
not here. This can be interpreted in the following way: 
At the transition, the spins that give rise to the edge 
states are no longer localized at the end of the chains. If 
they were free to move along the chain, they would give 
rise to a band of low-lying excitations and to a gapless 
spectrum. However, they move on a ladder, not on a 
chain, and this peculiar geometry opens a gap in their 
spectrum. 

For larger spins, a complete picture of the edge states 
cannot be obtained on the basis of exact diagonalizations 
because of the size limitation, and it will be necessary to 
use more sophisticated numerical approaches, which goes 
beyond the scope of the present paper. Still, as shown in 
the Supplemental Material, clear indications of changes 
in the edge states at the transitions could be obtained for 
spin-1, where, in the strong-rung limit, we expect to have 
two orphan spins 1 at each edge, hence 9 low-lying states 
below the gap (a singlet, a triplet and a quintuplet). 

As an alternative check of the valence-bond singlet 
picture, we have looked at the isotropic spin-3/2 ladder 
{J± — J\\)- 9 = 7r/4 lies in a phase with Berry phase tt, 
and the valence-bond singlet picture puts it in the phase 
with one singlet per rung. Since a spin 3/2 is able to build 
three valence-bond singlets, and since the coordinance of 
a 2-leg ladder is equal to 3, this phase is consistent with 
the valence-bond solid phase in which each bond (rungs 
and legs) carries one valence-bond singlet, i.e. in which 
the Berry phase of all bonds (and not just on the rungs) 
is equal to tt. So we have calculated the Berry phase 
on the leg bonds, and indeed it is equal to tt in that 
phase. So, it seems that the isotropic spin-3/2 ladder is 
a valence-bond solid. 

To further confirm this conclusion, we have studied the 
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Figure 3: (Color online) Top panels: Low energy spectrum of 
the spin-S =1/2 ladder with two types of boundary condi- 
tions (PBC and OBC with diagonal edges) for two values of 
the rung coupling: (a) in the rung phase {9 = 0.3 tt), where 
an edge state is clearly present for OBC with diagonal edges 
and (b) in phase without valence-bond singlets on the rungs 
(6* = O.Itt) where there is no edge state. Panel (c): Local 
magnetization in the sector 5*^1 = 1 for the spin S = 1/2 
ladder with diagonal edges. Inset: sketch of a cluster with di- 
agonal edges. Squares (circles) have been used for sites of the 
top (bottom) leg. (d): Correlation length as a function of 6 
as deduced from the exponential scaling of the edge state gap 
(4^) with the number of rungs Nr = N/2, and from fitting 
the magnetization profile (5*^). 

Hamiltonian H{a) = (1 - a)H^'"^'^''' + aH^^^^ that in- 
terpolates between the Heisenberg model and the AKLT 

modelE25l: 



^AKLT _ 


^11 XI XI P{i,a)Ai+i,o,) +J±^ 




i Q = l,2 


where Pi^j 


_ 11 1 27 a a 1 29 /o o \2 
^ 128 ' 160 ^« ^J ' 360 ^^« ^31 



E^( 



(^l),(^2), 



gp(Si • Sj) 



is the projector on total spin 3. When a = 1, the AKLT 
state is the ground state for any 9 in the range Q < 9 < 
7r/2 since Jy and J_l are strictly positive. As can be seen 
in Fig.|4J this AKLT phase at a = 1 smoothly evolves into 
the intermediate phase with a rung Berry phase tt and 
one valence-bond singlet per rung (0.203 t: < 9 < 0.308 tt) 
at a = 0. The leg Berry phase (not shown) is equal to tt 
in all this phase. This unambiguously demonstrates that 
the S=3/2 isotropic ladder is in the same phase as the 
AKLT model, hence that it is in the valence-bond solid 
phase. 

In the same spirit, the valence-bond singlet picture pre- 
dicts that the phase of the spin-2 ladder with two valence- 
bond singlets per rung should have one valence-bond sin- 
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Figure 4: (Color online) Rung Berry phase of the model that 
interpolates between the Heisenberg and AKLT models as 
a function of a for 8 sites. The AKLT valence-bond solid 
at a = 1 (sketched as an inset) evolves smoothly into the 
intermediate phase that contains the isotropic ladder at a = 
0. 



glet per leg bond, hence a Berry phase tt, a prediction 
confirmed by our numerical results. 

To summarize, we have shown that spin-S ladders with 
antiferromagnetic rungs undergo a cascade of topological 
phase transitions upon changing the ratio of leg to rung 
coupling. At the transitions, a topological invariant, the 
Berry phase on a rung, changes from to vr, and edge 
states appear or disappear, but the bulk gap does not 
close. The absence of gap closing might be taken as an 
indication that there is no true phase transition. How- 
ever, as first pointed out by Hatsugai|13]. the quantiza- 
tion of the Berry phase as equal to or tt module 2tt 
is a consequence of time reversal symmetry, and phases 
with Berry phase and tt are topologically distinct, the 
difference being protected by this symmetry in the spirit 
of Ref. 6. In that respect, the phases with the same 
Berry phase modulo 2 7r, i.e. with the same parity of the 
number of singlets in the valence-bond singlet interpre- 
tation, should not be considered as topological distinct 
phases, although the classification according to the num- 
ber of valence-bond singlets (and not simply its parity) 
is useful in discussing edge states. 

So the present results lead to the conclusion that, at 
the transition between different topological phases, the 
bulk gap does not necessarily close, even if edge states 
appear. What happens is that edge states indeed delo- 
calize at the transition, but the spectrum remains gapped 
because of another interaction, here the interaction be- 
tween the chains due to the ladder geometry. It will be 
interesting to see if similar topological phase transitions 
without gap closing appear in other contexts. 

We are indebted to E. Orignac and K. Penc for en- 
lightening comments on the fate of the edge states at the 
transition. This work has been supported by the Swiss 
National Fund and by MaNEP. 
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SUPPLEMENTAL MATERIAL FOR "TOPOLOGICAL PHASE TRANSITIONS IN 

SPIN-S LADDERS" 
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NUMERICAL CALCULATION OF THE BERRY 
PHASE 

In the main text, we have given a definition of the 
Berry phase as an integral of the Berry connection over 
a closed path: 



ij = (b A{if))d(l), 



(3) 



where A{(t>) == {GS{(t))\d^\GS{(j))). To calculate the Berry 
phase numerically, it is expedient, following Ref. 1, to 
introduce a gauge-invariant Berry phase on the lattice 
by discretizing in cf): 



K-l 



IK 



^arg(G^(0fc)|G5((/)fc+i)), 



k=a 



2Trk 



(4) 



The Berry phase is then given by 7 = limx->oo Jk- In 
practice, the convergence is achieved very fast. Through- 
out, we have used K — 16. 

Now, for large system sizes, it becomes numerically 
quite expensive to compute the Berry phase because one 
needs to diagonalize the system an important number of 
times (here 16), with just a change on one bond. How- 
ever, a change of Berry phase is associated to a level 
crossing, and this level crossing usually appears for sim- 
ple values of the twist parameter, = or = tt. Since 
we know that the gap never closes for = 0, the level 
crossing has to happen for (f> = tt, i.e. for anti-periodic 
boundary conditions. In other words, it is sufficient to 
compute the ground state energy of the system with anti- 
periodic boundary conditions, hence for a single value of 
0, and to determine the parameter at which a level cross- 
ing occurs. This analysis is done in Fig.[5]for A^ = 30 sites 
around the first transition at 9 « 0.26 tt. The systematic 
investigation of large system sizes has been performed 
with this method, and not by computing explicitly the 
Berry phase. 



FINITE-SIZE ANALYSIS OF THE TRANSITION 
POINTS FOR SPIN S = 1 AND S = 3/2 LADDERS 

In Fig. |6) we show the finite-size analysis of the transi- 
tions for spin S = 1 and S = 3/2 ladders. The results for 
5' = 1/2 can be found in the main text. The extrapolated 
values have been obtained by fitting the results with first 
or second order polynomials in 1/N. 

For spin S* = 1 , we did not observe any change of Berry 
phase for A^ = 6 and N ^ 10 for the first two transitions. 
This is an extreme case of the even-odd effect, which is 
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Figure 5: (Color online) Closing of the gap at the transition 
around 6 ~ 0.26 tt for the spin- 1/2 ladder with A'^ = 30 sites 
and twisted periodic conditions. A very precise estimate of 
the transition point can be simply obtained by a linear fit. 



generally present when the leg coupling is antiferromag- 
netic. So the fit has been performed using only the results 
for systems with an even number of rungs, i.e. N = 8, 
iV = 12 and AT = 16. 

For S = 3/2, the same is true for the transitions around 
e w 0.308 TT and 9 w 0.389 tt. For 9 « 0.203 tt, we clearly 
see a strong even-odd effect. As in the case of the spin- 
1/2 ladder, extrapolating the results separately for sys- 
tems with an even and odd number of rungs respectively 
leads to results in excellent agreement, with a difference 
of the order of 0.00005 in the present case. 



EDGE STATES FOR THE SPIN S = 1 LADDER. 

In this section, we summarize the results we have ob- 
tained regarding the appearance of edge states for spin-1 
ladders under open boundary conditions with diagonal 
edges. 

In the strong antiferromagnetic rung coupling limit, 
two isolated spins S = 1 appear at the edges, forming 9 
low lying energy levels, one quintuplet, one triplet and 
one singlet. Since we work in the subspace 5"*°' = 0, we 
expect to see only three states, namely the singlet, one of 
the triplets, and one of the quintuplets. The first one is 
the ground state. So, since we measure the energies with 
respect to the ground state energy, we expect to find two 
low- lying states. This is clearly the case in Fig. [t] (a), 
which corresponds to 9 — 0.45 tt, well inside the strong 
rung singlet phase with two valence-bond singlets per 
rung. 

The situation is quite different in the intermediate 
phase {9 = 0.3 tt) shown Fig. [71(b). Indeed, open bound- 
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Figure 6: (Color online) Finite-size analysis of the topological 
transitions in the spin-1 (left panel) and spin-3/2 (right panel) 
ladders. The fits have been performed with first and second 
order polynomials in 1/A'^ 



ary conditions with diagonal edges give rise to only one 
state below the gap. In addition, its energy is not so 
small. With the system sizes available, it is not possible 
to decide whether its energy goes to zero in the thermo- 
dynamic limit. To get additional information on this mid 
gap state, we have plotted in Fig. [t] (d) the local magne- 
tization in the sector S^^^. — 1. The local magnetization 
does not seem to decay exponentially. So it seems that 
the mid-gap state is not a well localized edge state. 

Finally, the situation changes again for weakly coupled 
chains, in which case open boundary conditions with di- 
agonal edges give rive to several mid-gap states (see Fig.lT] 
(c)). This can be understood as follows: In the limit of 
completely decoupled chains, there is an effective spin 
1/2 at each end, i.e. 4 in total. So we expect to find 
16 low-lying states in total. However, as soon as a small 
coupling between the chains is introduced, the two free 
spins 5* = 1/2 at one end of the ladder are expected to 
couple, leading to a finite excitation gap. This is con- 
firmed by the magnetization, which is no longer localized 
at the edges of the chain (see Fig. [t] (d)) presumably 
because the spin- 1/2 edges of the decoupled chains are 
bound into singlets. 



CORRELATION LENGTH FROM LOCAL 
MAGNETIZATION 

For the spin-1/2 ladder, when the system is in the 
strong rung coupling phase with one valence-bond sin- 
glet per rung, the correlation length can be obtained by 
fitting the decay of the lowest excitation energy with the 
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Figure 7: (Color online) Low-energy spectra of the spin 5=1 
ladder for different boundary conditions. In (a), we are in the 
rung singlet phase and we clearly see an edge states for OBC, 
diagonal edge. The fits use exponentially decaying functions 
in N. In (b) we are in the intermediate phase, (c) corresponds 
to two weakly coupled spin-1 chains. Valence-bond singlet 
sketches are attached to spectra (a) and (c). Dots denote 
spin-1/2, bold lines denote singlets, a single arrow indicates 
an isolated spin-1/2, and a double arrow stand for a spin-1. 
(d) Local magnetization in the sector St^t ~ 1- Squares and 
circles distinguish spins from top and bottom chains. 



size (see main text), or by fitting the decay of the local 
magnetization, as discussed in detail in this section. The 
local magnetization is defined as the expectation value of 
Sf in the ground state of the sector S^^^ = 1. The distri- 
bution of local magnetization is symmetric with respect 
to the vertical cut in the middle. Therefore only half the 
ladder can be taken into account for the fit. Moreover, 
in order to minimize the effect of the overlap of the left 
and right sides of the ladder, we do not consider the two 
rungs in the middle. Finally, we also did not include the 
magnetization of the first site to avoid the boundary ef- 
fect. Then a log-linear plot of the local magnetization 
leads to a linear curve with a very good accuracy, from 
which a correlation length can easily be extracted (See 
FigM. 
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Figure 8: (Color online) Log-linear plot of the local magne- 
tization for the 30-site spin-1/2 ladder with open boundary 
conditions and diagonal edges calculated in the ground state 
of the sector Stot = 1- The correlation length is extracted 
from a linear fit. 



PHASE TRANSITION AT 9 ^ iv AND 37r/2 

In the phase diagrams of the main text, we have indi- 
cated that the ferromagnetic phase extends from 6 — -k 
an 9 = 37r/2. That the gromid state is ferromagnetic 
in this range is clear since both couplings are ferromag- 
netic. To prove that the ferromagnetic phase does not 
extend beyond this range, we have compared the ground 
state energies in all sectors of S'fot' which are expected 



to be completely degenerate only if the ground state is 
ferromagnetic (see Fig. lol). 







Figure 9; (Color online) Ground state energies of a 16-site 
spin-1/2 ladder in the different sectors of Stot- The ground 
state is completely degenerate, hence ferromagnetic, for 9 G 

{7T,37v/2). 
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